(N 

O 
(N 

o 

Q 

m 



Theoretical Treatments of the Bound-Free Contribution and Experimental 
Best Practice in X-ray Thomson Scattering from Warm Dense Matter 
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Department of Physics, University of Washington, Seattle, WA 98195-1560 
(Dated: 14 December 2012) 

By comparison with high-resolution synchrotron x-ray experimental results, we assess several theoretical 
treatments for the bound- free (core-electron) contribution to x-ray Thomson scattering (i.e., also known as 
nonresonant inelastic x-ray scattering). Wc identify an often overlooked source of systematic error in the plane- 
wave form factor approximation (PWFFA) used in the inference of temperature, ionization state, and free 
electron density in some laser-driven compression studies of warm dense matter. This error is due to a direct 
violation of energy conservation in the PWFFA. We propose an improved practice for the bound-free term that 
will be particularly relevant for XRTS experiments performed with somewhat improved energy resolution at 
the National Ignition Facility or the Linac Coherent Light Source. Our results raise important questions about 
the accuracy of state variable determination in XRTS studies, given that the limited information content in 
low-resolution XRTS spectra docs not strongly constrain the models of electronic structure being used to fit 
the spectra. 
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I. INTRODUCTION 

The field of warm dense matter (WDM) rests in the 
transitional regime between traditional condensed phase 
systems and strongly-correlated, fully-ionized plasmas. 
As such, it draws from the complexity of both fields 
while showing its own special fundamental and, as wc 
show here, pragmatic challenges. A significant and 
growing literature exists on the electronic structure?i"— 
thermodynamics^ and hydrodynamics^ of WDM, in ad- 
dition to a range of applications in fusion energy science^ 
and laboratory astrophysics i^*^ Here, however, we inves- 
tigate a particular difficulty of the WDM regime: assess- 
ing the accuracy of the experimental determination of the 
basic state variables of the system, such as temperature, 
density and ionization state. Reliable inference of these 
quantities is central to the clearly-needed improvements 
in the equation of state of materials under, for exam- 
ple, the entire range of conditions leading from ambient 
matter to inertial confinement fusionJ^ 

The high optical opacity of WDM requires the use 
of penetrating probe radiation, i.e., x-ray photons with 
energies of a few to a few tens of keV. Unfortunately, 
with only few counterexamples^^ the inference of state- 
variables using these methods is limited by the degree 
of understanding of the electronic structure of WDM 
and its relationship to the state variables themselves. 
Faintly circular co- dependencies of this type are not un- 
common in emergent fields of experimental science (e.g., 
consider the many years of effort needed to establish ac- 
curate and precise pressure sensing in the Mbar range in 
opposed- anvil pressure cells^"— ), and a firm foundation 
for such methodologies can follow from any of several 
developments. Foremost among such developments are: 
experimental data of sufficient information content to it- 
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self strongly constrain the constituent theories for elec- 
tronic structure; broad programs to assess accuracy by 
cross-comparison of different metrologies; and, finally, an 
eventual comparison to international standards. The ex- 
perimental determination of state variables in the WDM 
regime is seeing only the earliest such examples, includ- 
ing notably the rare experiments using detailed balance 
in x-ray Thomson scattcringi^ or the recent checks in 
consistency between conclusions drawn from the clastic 
and inelastic components of x-ray scatteringJi 

The present paper reports a first step in evaluating the 
accuracy, rather than precision, of the methods used for 
state variable determination in the WDM regime. To 
this end, we investigate the various available treatments 
for the core-electron, or bound-free, contribution to x-ray 
Thomson scattering (XRTS, also called nonresonant in- 
elastic x-ray scattering, NIXS^^) with a special emphasis 
on the plane-wave form-factor approximation (PWFFA) 
of Schumacher, et alJ^. We will show that this approxi- 
mation, which has seen extensive use in XRTS studies of 
shock-compressed matter ;^ ''i^"" — is fundamentally flawed 
and presents a source of systematic uncertainty in in- 
ferred quantities. 

Based on these observations wc come to three main 
conclusions. First, looking to the near future, when 
XRTS studies of WDM with improved energy resolution 
will be performed at the Linac Coherent Light Source^^ 
and the National Ignition Facility^, the errors implicit 
in the use of the PWFFA must be avoided if physically 
meaningful information on the equation of state in the 
WDM regime is to be determined. Second, while it is 
important to note that a faulty theoretical treatment 
has been used, and that some rcevaluation of experi- 
mental results may be called for, the more lasting con- 
clusion is that the information content in the measured 
XRTS spectra for WDM has been insufficient to alert 
the experimenters to the presence of an unphysical model 
for the electronic structure. This strongly suggests the 
need for cross-comparison with alternative methods of 
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WDM state variable determination, e.g., x-ray fluores- 
cence thermometrj*22r— . Third, we find that stronger 
connections between the synchrotron x-ray and WDM- 
XRTS communities provide important experimental and 
theoretical synergies. The wealth of very high-resolution 
studies at synchrotron light sources of both the free-free 
(valence)^"— and bound-free (core) ^^i^^" — contributions 
to XRTS provide important benchmarks both for com- 
parison to WDM-specific theory and also for validation of 
experimental protocol including, e.g., instrument-specific 
backgrounds. Further, as we have illustrated here, there 
will be cases where theoretical methods already in use 
for synchrotron studies may be beneficially transported 
to WDM-XRTS studies. 

In Sect. In] we discuss four theoretical treatments 
of bound- free XRTS: the impulse approximation (lA), 
which is valid at large energy transfer; a hydrogcnic 
model (HM) ; an extension of the lA to incorporate bind- 
ing energies (PWFFA); and a real-space Green's function 
method (RSGF). The first three are essentially atomic 
(with varying degrees of approximation), while the lat- 
ter treats the condensed solid, and is based on methods 
broadly used for many years in the interpretation of sev- 
eral x-ray spectroscopic techniques 

Since high-resolution XRTS data from WDM at known 
thermodynamic conditions is currently unavailable, we 
instead compare each of the above theories with very 
high-quality XRTS spectra collected under ambient con- 
ditions at a synchrotron x-ray source. This provides a 
baseline validation for the core contribution: a theoreti- 
cal treatment which fails under these conditions is certain 
to form a weak foundation when including the further 
complexities of continuum lowering and partial ioniza- 
tion present in WDM. 

Experimental details are described in Sect. IIIII and 
the comparison of theory and experiment is made in 
Sect. IIV Al The relative success of even atomic treat- 
ments at describing the condensed solid suggests that 
these methods should be extensible to the WDM regime 
with only minor modifications. However, we find that 
the PWFFA, which has been used for a few years in the 
interpretation of WDM measurements jiii^"— is in stark 
disagreement with the ambient experimental data. In 
Sect. IIVBI we show that this disagreement is due to in- 
ternal inconsistency in the PWFFA that leads to unphys- 
ical results. Next, in Sect. IIV Cl we consider the implica- 
tions of using the PWFFA to model the bound-free con- 
tribution to WDM XRTS data; namely, the likelihood 
of previously-undiagnoscd systematic errors in extracted 
thermodynamic quantities. These observations then mo- 
tivate a discussion of best future practice in Sect. IIVDI 
after which we conclude in Sect. El 



II. THEORY 

The fundamental observable in XRTS/NIXS is the dy- 
namic structure factor S{q, w) , which separates into inde- 



pendent contributions from electrons in different shells. 
We will focus on the contribution from tightly bound 
core electrons, i.e., the bound-free contribution. The the- 
oretical description of bound-free XRTS begins with the 
Kramers-Heisenberg formula for the first-Born approxi- 
mation to the double-differential scattering cross-section 
(DDSCS)>i^ 
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Here, a;i_2 and ei^2 are initial and final photon energies 
and polarizations; rg is the Thomson scattering length; 
j/),!-?") are initial and final many-body states with ener- 
gies Ej.Ep] Pi[T) is the temperature-dependent Boltz- 
mann factor; q is the momentum transfer; w = wi — a;2 
is the energy transfer and j indexes individual electrons. 
In this and subsequent formulae, we use Hartree atomic 
units {Ti = me = 1). 

In the independent particle approximation, Eq. ([2]) can 
written as 

S{q,oj) = ^n,Si{q',uj) 

i 

sM ^) - E(i - I (/I I*) r - - 
/ 

where are initial and final state single particle 

states with energies Ei,Ef and thermally-averaged oc- 
cupation numbers n^, nj. 



A. The Impulse approximation 

In the limit of large energy-transfer uj relative to the 
initial state binding energy Eb, known as the impulse 
approximation (lA), the XRTS spectrum is completely 
determined by the initial-state electronic momentum dis- 
tribution; the binding energy of the scattering electron 
plays no rolci^ 

For Lj ^ Eb, only unoccupied final states contribution 
to Eq. (I3|), so we set n/ = 0. Next, following Eisen- 
berger and Platzman^^ we expand the (5-function using 
the standard Fourier representation 



(4) 



After rearranging slightly and using the fact that \i) and 
I/) are cigcnstatcs of the single-particle Hamiltonian H, 
we find 



SM.^) = j ^E 



dt 
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In the second line, we have used completeness to remove 
the sum over final states. The I A corresponds to replac- 
ing H by the free-particle Hamiltonian Hq, which can be 
justified in the limit of {oj/EbY ^ 1 (see Section III of 
Ref. [56| ). After inserting a complete set of momentum 
eigenstates and integrating over the direction of momen- 
tum, we obtain 



S^{q,uj) = (27r/q) 



w/q~q/2\ 



pdppiip), 



(6) 



where pi{p) = (27r)~'^| {i\p) p is the initial-state momen- 
tum density, which is here assumed to be isotropic (e.g. 
s-shell, or sum over a filled subshell). 

This formula can be interpreted as describing XRTS 
from a gas of free electrons with the same initial-state 
momentum distribution. The scattering spectrum con- 
sists of a line centered at the free-particle Compton shift 
ujc ~ 9^/2 that is Doppler broadened by the projection of 
the momentum distribution along the direction of q. The 
integral over the momentum distribution in Eq. (|6]) sim- 
ply counts electrons with a given momentum-projection 
Pq determined by the energy transfer. 

We now turn to methods that take the binding energy 
into account. 



B. Hydrogenic Model 

It is possibly to analytically evaluate Eq. ([3]) using 
hydrogenic initial and final states. For a Is initial state, 
the result is^ 



Siiq,uj) 



with 
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X [fcV + (l/3)fcV(l+p2a2)] 

X [(fc2a2 + l-pV)2 + Va2]-^ 
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Here, p is the final-state momentum, a = l/Z where Z 
is the effective nuclear charge^, and Eb = 1/2Z^ is the 
binding energy. In this expression, contributions from 
bound final states have been neglected. Expressions for 
shells other than the Is are included in Schumacher, et 
aliiS , where this method is referred to as the hydrogenic 
form-factor approximation. We will, however, refer to 
this simply as the hydrogenic model (HM). 



C. Plane-wave form-factor approximation 

The plane-wave form-factor approximation (PWFFA) 
is an attempt to improve the lA by including the bind- 



ing energy in the kinematics. The final states are as- 
sumed to be momentum eigenstates, which is concep- 
tually appealing in the context of dense plasmas where 
the jellium model has found wide application. However, 
as we demonstrate in Sect. IIVBI a fundamental diffi- 
culty arises: this approximation effectively evaluates the 
initial-state energy using the atomic Hamiltonian, while 
evaluating the final-state energy using the free-particle 
Hamiltonian. This inconsistent treatment violates energy 
conservation, resulting in violations of the Bethe /-sum 
rule and deviations from experimental results that, al- 
though small in the original context of gamma-ray scat- 
tering, are quite large under the kinematic conditions 
typical of XRTS measurements. It is the use of the 
PWFFA in the interpretation of recent XRTS experi- 
ments on WDMiii^"— that motivates the present paper. 

Schumacher's derivation of the PWFFAii^ makes the 
following assumptions: 

I/) = \p + q) 

d^p 



E 



(27r)3 



Ef - - 

Ei = —Eb 



{P + q)^ 



(9) 



where Eb is the initial-state binding energy, p is the 
initial-state momentum and p + q is the final-state mo- 
mentum. Assuming that T = and applying ([9]) to ([3|), 
we find 



S,{q,uj) = 



^ I {p + q] e'^~-^ \^)\' 5{Ep+^ + Eb - u) 
d^p 



' S{Ep+^+ Eb - w) 



(27r)3 ' 

= J d'^pp,{p)5{Ep+^+EB-Lo) 



(10) 



In the second line we use the fact that e*'' '' is a momen- 
tum translation operator. The last line uses the defi- 
nition of the momentum density Pi(p) = {2'k)~^ \ {p\i)\ ■ 
Furthermore, if we again restrict ourselves to an isotropic 
momentum density (e.g., s-shell, or sum over a filled sub- 
shell), we can perform the angular integrals to obtain 



l.^2(L0-B)+q 

St{q,uj) = {2Tr/q) pdppi{p) (11) 

J\^2{L0-B)-q\ 

This expression, along with ([T]) differs from 
Schumacher'si^ Eqs. (5) and (21) only in that it does not 
contain the relativistic prefactor -^/l + (p/mc)-^, which 
for the experimental conditions under consideration 
differs negligibly from unity. 

Eq. ([TT|) has the same form as the lA expression 
Eq. (|6]), differing only by the bounds on the integration 
over the momentum density. Given this similarity and 
the well-tested validity of the lA in its regime of appli- 
cability, one would expect that for to » Eb the PWFFA 
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would reproduce the lA. This is, however, not the case 

— despite claims in the literature to the contraryJ^i^ 
We will return to this point in Sect. lIV| after comparison 
with experiment. 

D. The real-space Green's function method 

The prior methods have all treated XRTS for an iso- 
lated atom to various degrees of approximation. The 
next method we describe treats an arbitrary cluster of 
atoms using a real-space Green's function (RSGF) for- 
malism implemented in recent versions of the x-ray spec- 
troscopy code FEFF.— This formalism, which can treat 
complex, aperiodic systems, has been extensively applied 
to condensed-matter systems, where XRTS/NIXS pro- 
vides a bulk-sensitive alternative to, and extension of, 
soft x-ray absorption spectroscopy , ^^^ ^ The RSGF ap- 
proach has recently been extended to treat the valence 
contribution)^ 

Starting with a description of atomic species and lo- 
cations, an effective one-particle Green's function for the 
valence electrons in the cluster of atoms is calculated in 
the muffin-tin approximation, including the effects of full 
multiple scattering^. This Green's function implicitly 
contains the excited electronic states that are the final 
states in the scattering experiment. In terms of a spectral 
density matrix defined by p{E) — J2f I/) (/I ~ Ef), 
which is related to the Green's function by {f] p{E) \ f') = 

— (I/tt) ImG'(r', f, -E), Eq. ([3|) can be recast as 

S,{q,u:) = (*| e-''!-^^'Pp{E)Pe'^-^^\i) . (12) 

Here E = lo + Ei \s the photoclcctron energy and P 
projects the final states (which are calculated in the 
presence of a core hole) onto the unoccupied states of 
the initial-state Hamiltonian (which has no core hole)»^ 
The Green's function can be separated into contributions 
from the central atom and from scattering off other atoms 
in the cluster. Likewise, the dynamic structure factor can 
be factored as 

SM^)^So{q.Lo)[l + Xq{^% (13) 

where So{q,uj) is a smoothly varying, isotropic atomic 
background and Xq is the fine structure due to all orders 
of photoclcctron scattering from the environment 
Implicit in the fine structure is information about 
nearest-neighbor distances and thus also density^. How- 
ever, at the poor experimental resolution typical of WDM 
measurements^^, this structure will be washed out. Thus, 
for our purposes, we will include only the atomic back- 
ground contribution, SQ{q,io). 

III. EXPERIMENT 

Although we are ultimately interested in the elevated 
temperatures and densities of WDM, it is important 
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FIG. 1. (Color online.) XRTS from polycrystalline beryllium 
under ambient conditions at a fixed 171° scattering angle and 
9890-eV scattered photon energy^. Here, the energy transfer 
ijj is the difference between the incident and scattered photon 
energies. In the upper panel, the data are shown along with 
a combined real-space Green's function (RSGF) valence and 
core calculation. The data have been scaled as described in 
the text. In the lower panel, the valence contribution and 
linear background have been subtracted to give the core con- 
tribution alone. 



to first validate theoretical methodology against spec- 
tra taken under known thermodynamic conditions and 
at higher resolution than presently typical of WDM ex- 
periments. To this end, experimental XRTS data for 
polycrystalline Be at ambient temperature and pressure 
were collected using the lower energy resolution inelas- 
tic x-ray (LERIX) spectrometer at beamline 20-ID of 
the Advanced Photon Source.— Scattered photons with 
ijJ2 = (9891.7 ± 0.2) cV were analyzed by a single spher- 
ically bent Si crystal located at a fixed 171° scattering 
angle while scanning the incident photon energy. From 
the elastic peak width the total instrumental resolution 
was determined to be 1.3 eV. The data, which have pre- 
viously been reported^, are shown in Fig. [TJ The graph 
is labelled S{qg,uj) to indicate that the data are collected 
at fixed scattering angle, and thus g is a weak function 
of w. 

After normalizing to the incident flux, a small linear 
background and a single scale factor were fit in order 
to match the RSGF core calculation in the tail region 
(800 <u < 1500 eV) and thus put the data into absolute 
units. The results of this fit have been used to scale the 
experimental data in Fig. [TJ Additionally, theoretical 
core and valence calculations, the fit linear background, 
and the sum of these three are shown. 
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FIG. 2. (Color online.) Comparison of extracted core-shell 
XRTS with theoretical calculations using the RSGF, the 
hydrogenic model (HM), impulse approximation (lA), and 
plane- wave form-factor approximation (PWFFA) . The energy 
transfer lo is the difference between the incident and scattered 
photon energies. All calculations are in absolute units. Ver- 
tical guides are shown at the Is binding energy (112 eV) and 
the free-particle Compton shift (396 eV). 



The experimental generalized oscillator strength 
J{2/q'^)ujS{q,uj)duj matches that for the combined theo- 
retical RSGF spectrum to within 1%. Since the measure- 
ment was performed at fixed scattering angle, this value 
is slightly larger than the Bethe /-sum rulei^^iS value 
of N=4 (which only holds for experiments performed at 
fixed q). 

The theoretical valence profile, calculated from the 
RSGF^ was then subtracted to obtain the experimen- 
tal core profile, shown in the lower panel of Fig. [1] The 
small peak visible for 550 < a; < 650 eV is a result of 
the theoretical valence profile underestimating the actual 
contribution in this region, as was also seen in compar- 
isons with higher momentum-transfer data (see Figs. 4 
and 5 of Mattern, et al^^). 



IV. RESULTS AND DISCUSSION 

A. Comparison of Experiment and Theory 

The extracted experimental core profile is compared 
in Fig. [2] to to each of the four theoretical calculations 
discussed in Section [TTl All calculations are in absolute 
units and have taken the weak dependence of g on w into 
account. Vertical guides are included at the Is binding 
energy (112 eV) and the free-particle Compton shift (396 
eV). 
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FIG. 3. (Color online.) On the left are shown the integra- 
tion bounds for both the PWFFA (upper and lower bounds) 
and lA (lower only, upper bound is oo) calculations. On the 
right, plotted vertically, is the integrand pp{p). Both the off- 
set of the peak by the binding energy (112 eV) and the lack of 
convergence of the PWFFA to the lA at large lu are apparent. 



The lA and PWFFA calculations use the ground-state 
Dirac-Fock Be Is wavefunction calculated using FEFF's 
atomic solver as the initial state, and thus only differ in 
their treatment of the final states and energy conserva- 
tion. For the RSGF calculation, the Dirac-Fock wave- 
function is calculated in the presence of a Is core-hole. 
We have included only the atomic background contribu- 
tion So{q, Lu). The fine structure visible in the experimen- 
tal data for uj < 200 eV is not included, although it has 
been treated elsewherej^ The HM calculation uses hy- 
drogenic wavefunctions with an effective nuclear charge 
(Z=3.685)^ for the initial and final states. 

We focus on three regions for comparison: the vicini- 
ties of the X-edge binding energy (w ^ 112 eV), the peak 
(w - 396 eV), and the tail (w > 600 eV). The RSGF cal- 
culation matches the data reasonably well in all three 
regions (with the exception of immediately above the K 
edge, where interference effects have been omitted). The 
HM accurately describes the peak and tail regions, but 
shows a large deficit above the experimental binding en- 
ergy due to the larger binding energy of the hydrogenic 
state. The lA, which ignores the binding energy, has an 
unphysical tail at low energy transfers. The peak region 
is reasonably well described by the I A, while the high-w 
tail region is quite accurate. This is expected since the 
conditions of applicability of the lA are well satisfied for 
large energy transfer. 

By contrast, although the PWFFA of Schumacher, et 
alJ^ vanishes below K edge, it exhibits only a gradual 
onset and further shows strong quantitative and qualita- 
tive disagreement with the experimental data everywhere 
else. Given its application in the interpretation of sev- 
eral XRTS experimentsiii^"— on WDM, and its evident 
failure to describe high-resolution synchrotron measure- 
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mcnts, we now turn our focus to the PWFFA. We will 
first look in more detail at the source of the approxima- 
tion's error, and then briefly discuss the possible impli- 
cations for interpretation of experiment and future best 
practice. 



B. A closer look at the PWFFA 

Although others have previously observed that the 
PWFFA gives results with unphysical features that are 
in disagreement with experimental data,— >^ we are un- 
aware of a discussion of the origin of the approximation's 
inconsistency. We now consider this point in detail. 

An alternative route to obtain the PWFFA is to follow 
the I A derivation up to Eq. ([5|). At this point, if one 
makes the ad hoc approximation of replacing only the 
second H hy Hq, 



dt 
2^ 



(14) 



then, instead, the PWFFA result ([TT]) follows. Thus, the 
assumptions Q correspond to making the uncontrolled 
approximation of evaluating the initial-state energy using 
H and the final-state energy using Hq. This effectively 
violates energy conservation, opening the possibility of 
unphysical results. 

As we mentioned in Sect. Ill CI the lA ^ and 
PWFFA (fTTjl results differ only in the bounds of the in- 
tegration over the momentum density. In left panel of 
Fig. [31 we show the integration bounds as a function of 
u) for both theories. For the lA, there is no upper bound, 
and for the PWFFA the upper bound is only relevant for 
small uj, beyond which it is well above the integrand's 
region of support. The right panel shows the integrand 
(rotated so that the abscissa runs vertically). The scat- 
tering spectra can be seen to peak when the lower inte- 
gration bound vanishes. For the PWFFA, this is offset 
to higher energy transfer by the binding energy of the 
initial state. Furthermore, the failure of the PWFFA to 
reduce to the lA at large uj can be clearly seen here. The 
offset of the PWFFA peak relative to the lA appears to 
be in conflict with the calculations presented in Fig. 3 of 
Riley, et alM. 

Given the much lower energy resolution typical of 
WDM experiments, it is natural to ask whether the dis- 
crepancies seen above are relevant in that context. In 
Fig. m we compare RSGF, lA and PWFFA calculations 
at the representative experimental conditions of Fort- 
mann, et alJ^. The lA calculation has been truncated at 
the empirical binding energy. In addition to the absolute- 
unit PWFFA calculation, we have included a curve scaled 
to match the /-sum of the other theoretical curves. The 
unbroadened calculations are shown in the upper panel. 
The curves in the middle panel are broadened by 115- 
eV (FWHM) to match the experimental resolution of 
Ref.[l3. Finally, for illustrative purposed, the lower panel 
contains curves broadened by 500 eV. Even in the latter. 
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FIG. 4. (Color online.) Theoretical core profiles for the ex- 
perimental conditions of Fortmann, et al.— . The inaccuracy 
of the PWFFA remains, even after substantial broadening. 
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FIG. 5. (Color online.) Theoretical (RSGF) XRTS from 
polycrystalline Be compressed to 3.6 x ambient density. The 
valence contribution is broader than under ambient condi- 
tions (c/. Fig. [Tl top panel). Subsequently, the core contribu- 
tion (assumed here to be independent of density) is relatively 
larger in the region of the peak. 



admittedly extreme case, the PWFFA oS'set and sub- 
sequent overestimation of the high-w tail is still quite 
prominent. The ad hoc /-sum scaling results in a tail 
that is only slightly overestimated, but at the cost of an 
extreme deficit beneath the peak. 

At the higher densities typical of WDM, the core con- 
tribution becomes even more important. As the density 
is increased, the Fermi level increases relative to the bot- 
tom of the valence band resulting in a broader valence 
contribution. The core wavefunction, on the other hand, 
is only weakly dependent on density (at least for mod- 
est compression, where the cores from neighboring sites 
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have negligible overlap). Subsequently, as shown in Fig. [5] 
(c/. Fig.[TJ top panel), the core contribution is relatively 
larger in the peak region, increasing the importance of a 
numerically accurate theoretical treatment. 



C. Implications 

We now turn to implications of an incorrect core treat- 
ment on the interpretation of XRTS spectra from WDM. 
It is important to recognize the difficulty of these exper- 
iments and their analysis. The intrinsic width of back- 
lighter x-ray sources fundamentally limits the energy res- 
olution obtainable in this measurement technique. The 
low flux and need for single-shot measurement requires 
the use of low-resolution spectrometers with limited spec- 
tral range further decreasing the resolution while also 
complicating background characterization. This uncer- 
tainty in the background subtraction makes /-sum nor- 
malization especially difficult, and thus the spectra often 
can not be reliably placed into absolute units. Further- 
more, the highest likelihood background is necessarily 
dependent upon assumptions made about the core con- 
tribution to the spectrum. 

The complicated interplay between the various degrees 
of freedom present in such fits makes it difficult to state 
the exact implications of using the PWFFA in the ex- 
traction of thermodynamic state variables in published 
work.— i22r— It is likely that, in order for a good fit in 
the high-w tail to be obtained, the ionization state must 
be overestimated. This could explain the discrepancy 
noted by Fortmann, et al.^^ between their best-fit ion- 
ization state found using the PWFFA and earlier work 
that appears to use the HMi2£i21 Beyond that, the net ef- 
fect on extracted thermodynamic parameters is unclear. 
However, due to this previously undiagnosed system- 
atic uncertainty, re-evaluation of existing experimental 
datapiii^Sr— using a more appropriate core calculation and 
a maximum likelihood treatment of the background is 
necessary. 

D. Future Practice 

The limitations of the PWFFA bring up the question 
of best future practice for fitting the core-shell XRTS 
from WDM. This will become particularly important 
when higher resolution WDM- XRTS experiments are 
performed at the Linac Coherent Light Source and the 
National Ignition Facility. If such experiments are to 
reach their full scientific potential, errors on the scale 
of those given by the PWFFA must be avoided. An 
ideal treatment would include self-consistent determina- 
tion of occupied and unoccupied electronic states includ- 
ing condensed-phase effects. Also, the decrease of ion- 
ization potentials with increased density (i.e., continuum 
lowering) should be either implicitly present in the calcu- 
lation, or tunable using models from plasma physicsi^"— 
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FIG. 6. (Color online.) Comparison of theoretical core contri- 
bution to Be K-edge XRTS calculated using RSGF and trun- 
cated lA methods as a function of momentum transfer. Ver- 
tical arrows are located at the free-particle Compton shifts. 

Of the methods we have presented, the RSGF calcu- 
lation most closely describes the ambient experimental 
data. Condensed-phase effects are included explicitly in 
final states. Although a frozen atomic core wavefunc- 
tion is used, this is a common feature of all techniques 
under consideration, and should be sufficient at modest 
densities where core overlap is expected to be negligible. 
The primary limitation of the RSGF approach is that it is 
currently unclear how to incorporate continuum-lowering 
effects. 

Alternatively, one can use an ad hoc modification of 
the low-energy-transfer tail of the lA 

Su-lA{q,uj) = S{q,u) (^1 - -^^—Ly—^ . (15) 

We will refer to this approach, which for T = sim- 
ply truncates the spectrum below binding energy, as the 
truncated I A (tr-IA) . This allows straightforward applica- 
tion of continuum-lowering models to adjust the binding 
energy. A similar approach, using screened hydrogenic 
wavefunctions for the initial state instead of Dirac-Fock 
is discussed in Gregori, et ali^, but only in the context of 
smaller where the core contribution is relatively small. 
In Figs. [BJ [71 and [51 we explore the accuracy and appli- 
cability of the tr-IA at T=0 by comparing it with RSGF 
calculations for a range of momentum transfers for the K 
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FIG. 7. (Color online.) Same as Fig.[6l except for Al Li- and 
L2,3-edg6 XRTS. 



the truncation of the low uj tail results in slight /-sum 
violations (< 5% for Al, q>& A~^). We also note that 
for the Al L shell, the tr-IA and RSGF calculations differ 
for the individual subshells (Fig. |S1 upper two panels). 
However, agreement is recovered after combining to form 
the total L-shell contribution (Fig. [5] lower panel). 

Recently, another approach to modeling XRTS from 
WDM has been discussed by Johnson, et They use 
an average- atom model, which gives a significant im- 
provement in the treatment of the free-electrons com- 
pared to a simple jellium model. Unfortunately, the 
average-atom binding energies disagree significantly with 
experiment, limiting the accuracy of the bound-free con- 
tribution to the XRTS spectrum. Johnson, et ali2^ also 
consider applying the PWFFA to the average-atom Is 
state for Be and find similar qualitative discrepancies as 
we have discussed here. 

In summary, for XRTS experiments with modest en- 
ergy resolution at high momentum transfer it should be 
sufficient to treat the bound-free contribution with a 
truncated lA, where the truncation energy is adjusted to 
include the effects of continuum lowering. However, the 
I A ceases to be accurate at lower momentum transfers. 
If continuum-lowering shifts and temperatures are neg- 
ligible compared to the desired energy resolution, then 
the RSGF approach can be immediately applied at any 
momentum transfer. Further investigation is needed to 
determine if continuum lowering can be calculated or in- 
cluded empirically within the RSGF framework. 
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FIG. 8. (Color online.) Comparison of RSGF and lA cal- 
culations for Li (2s) and L2,3 (2p) subshells, along with the 
combined spectra at high q. While the RSGF and lA differ for 
individual subshells, agreement is recovered for the combined 
spectra. 

shell of Be and L shell of Al. Note that all calculations 
are in absolute units. 

As long as the free-particle Compton shift (shown by 
vertical arrows) is a few times the edge energy, the tr-IA 
is in reasonable agreement with the RSGF calculation. 
Since the lA satisfies the /-sum rule by construction. 



V. CONCLUSION 

We have discussed several techniques for calculating 
the core-shell contribution to XRTS and compared with 
experimental data collected from polycrystalline Be un- 
der ambient conditions. Of the techniques considered, 
the real-space Green's function method best describes 
the data. However, a simple ad hoc truncation of the 
impulse approximation is reasonably accurate at higher 
momentum transfers and allows more straightforward in- 
clusion of continuum lower effects. The accuracy of this 
truncated lA as a function of q has been explored by 
comparing with RSGF calculations for both the Be K- 
shell and Al L-shell. On the other hand, the plane- 
wave form-factor approximation, which has been used in 
the interpretation of several WDM experimentsiii^Sr— is 
quantitatively and qualitatively inaccurate due to an in- 
consistent treatment of the single-particle Hamiltonian. 
Re-evaluation of the experimental data using a more ac- 
curate core calculation and maximum likelihood back- 
ground subtraction is recommended. More importantly, 
an accurate treatment of the bound-free XRTS from 
WDM will be necessary when higher resolution experi- 
ments are performed at the Linac Coherent Light Source 
and the National Ignition Facility. We believe we have 
also motivated the need for cross-method comparisons 
and the usefulness of exchange of both theoretical and 
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experimental techniques between the condensed-matter 
and dense-plasma communities. 
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